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ABSTRACT 

We present parameter estimation forecasts for future 3D cosmic shear surveys for a 
class of Unified Dark Matter (UDM) models, where a single scalar field mimics both 
Dark Matter (DM) and Dark Energy (DE) . These models have the advantage that they 
can describe the dynamics of the Universe with a single matter component providing 
an explanation for structure formation and cosmic acceleration. A crucial feature of 
the class of UDM models we use in this work is characterised by a parameter, Coo (in 
units of the speed of light c = 1), that is the value of the sound speed at late times, 
and on which structure formation depends. We demonstrate that the properties of the 
DM-like behaviour of the scalar field can be estimated with very high precision with 
large-scale, fully 3D weak lensing surveys. We found that 3D weak lensing significantly 
constrains Coo, and we find minimal errors Acoo = 3.0 ■ 10~ 5 , for the fiducial value 
Cqo = 1-0 • 1CU 3 , and Acoo = 2.6 • 1CU 5 , for Coo = 1.2 ■ 1CU 2 . Moreover, we compute the 
Bayesian evidence for UDM models over the ACDM model as a function of Coq. For 
this purpose, we can consider the ACDM model as a UDM model with = 0. We find 
that the expected evidence clearly shows that the survey data would unquestionably 
favour UDM models over the ACDM model, for the values Coo > 1CU 3 . 

Key words: gravitation, gravitational lensing, cosmology: theory - observations - 
dark matter - dark energy - large-scale structure of Universe. 



"-H . 1 INTRODUCTION 



Provided that General Relativity itself is accurate, a number 
of cosmological observations, e.g. the dynamics of galaxies 
and galaxy clusters and the Large-Scale Structure (LSS), 
provide us with ample evidence that most of the matter in 
the Universe is not mad e up of familiar baryonic matter, but 
of Dark Matter (DM) JZwickvlll933l. Il937l; iDodelson et aT 



2001] : lHawkins et all 12003! : ISpergel et all 120071 ; iRiess et al 
20071 ). Moreover, the largest fraction of the energy budget 
in the present Universe is occupied by Dark Energy (DE), 
which seems to be accelerating the expansion of the Universe 
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(IRiess et al. lll9 98l;lKnop et al.ll2003l ; IRiess et al.ll2004 120071 ; 
lKomatsu~e t al. 2009). The current concordance model, with 
radiation, baryons, cold DM and DE (in the form of a cos- 
mological constant A) is known as th e ACDM mode l. 

It has been demonstrated (|Albrecht et al . 2006; 
IPeacock et al.l [§006) that the nature of the dark compo- 
nents of the Universe can be constrained to a high de- 
gree of accuracy by using wide and deep imaging sur- 
veys; weak lensing, in which the shear and redshift infor- 
mation of every galaxy is used, has the potential to con- 
strain the equation of state of such dark components by 
using surveys such as EuclicQ dRefregier fc Douspisl [2003 ; 
Refregier et al.ll2010h or Pan-STAR R£| | (|Kaiser et al.ll2002l ; 



Kaiser fc Pan-STARRS Team! [20021 ) . As a direct probe of 



1 http: / /sci.esa.int/science-e/www/area/index.cfm?farcaid=102 

2 http://pan-starrs.ifa.hawaii.edu 
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the mass distribution, gravitational lensing is an excellent 
tool for cosmological parameter estimation, complementing 
Cosmic Microwave Background (CMB) studies. One of the 
most useful manifestations of gravitational lensing by in- 
tervening matter is the alignment of nearby images on the 
sky. Detection of DM on large scales through such cosmic 
shear measurements - the small, coherent distortion of dis- 
tant galaxy images due to the large-scale distribution of mat- 
ter in the cosmos - has recently been shown to be feasible. 

At a statis tical level, it has been shown (|Hu fc Tegmarkl 
ll999l ; lHu1ll999h that there is some extra information on cos- 
mological parameters which can be gained by dividing the 
sample into several redshift bins; this technique is known as 
weak lensing tomography. However, a more comprehensive 
repr esentations of the shear field can b e called 3D weak lens- 
ing (lHeavensll2003l;ICastro et al.ll2005l ; iHeavens et aHl2006l ; 
iKitching et al. 20071) . in which, by using the formalism of 
spin- weighted spherical harmonics and spherical Bessel func- 
tions, one can relate the two-point statistics of the harmonic 
expansion coefficients of the weak lensing shear and conver- 
gence to the power spectrum of the matter density perturba- 
tions. Such a tool is relevant in view of the present and next 
generations of large-scale weak lensing surveys, which will 
provide distance information of the sources through photo- 
metric redshifts. 

Recently, rather than considering DM and DE as two 
distinct components, it has been suggested the alternative 
hypothesis that DM and DE are two states of the same fluid. 
This has been variously referred to as "Unified Dark Matter" 
or "quartessence" models. Compared with the standard DM 
plus DE models (e.g. ACDM) , these models have the advan- 
tage that we can describe the dynamics of the Universe with 
a single scalar field which triggers both the accelerated ex- 
pansion at late times and the LSS formation at earlier times. 
Specifically, for these models, we can use Lagrangians with 
a non-canonical kinetic term, namely a term which is an ar- 
bitrary function of the square of the time derivative of the 
scalar field, in the homogeneous and isotropic background. 

Originally this method was proposed to have 
inflation driven by kinetic energy, called fc-inflatio n 
l|Armendariz-Picon et al ]| 19991 ; IGarriga fc Mukhanovll 999). 
to explain early Universe's inflation at hi gh energies. 
Then this scenario was app li ed to DE JChiba et al.l 
|2000| ; Ide Putter fc Linderl 120071 ; iLinder fc Scherrerl 120091 ) 
In particular, th e analysis was extended to a more gen- 
eral Lagrangian l|Armendariz-Picon et ail I2OO0I, l200ll) and 
this scenario was called fc-essence (see also IChiba et al 



2000 1; iRendanl 120061 ; ILi et all l2006l: ICalcagni fc Liddk 



2006; Babichevl 120061; iFang et al. 20071; iBazeia et al 



Kangetall 120071 ; iBabichev et all 120081 ; iBabichevI 



2007 



2008 



Ahn et al.ll2009h . 



For zmodels, several adiabatic or, equivalently, purely 
kinetic models have been inves tigated in the literature : 
the generalised Chaplygin gas jKamenshchik et all l200ll ; 
Bilic et al.ll2002l ; lBento et al.|[2 002; Carturan & Fincllil l2003l ; 



Sandvik et al.ll2004 ) , the single dark perfect fluid with a sim- 



ple two-parameter barotropic e quation of state (|Balbi et alj 
l2007l ; lOuercellini et al. I l2007l ; iPietrobon etail 120081 ) and 



the purely kinetic models studied b y Scherrei l|2004h , 



iBertacca et~ai1 (|2007l ). IChimento et~aT1 (|2009l ). Alternative 
approaches have been proposed in mod els with can onical 
Lagrangians with a complex scalar field (|Arbevll2006l ). 



One of the main issues of these UDM models is to see 
whether the single dark fluid is able to cluster and produce 
the cosmic structures we observe in the Universe today. In 
fact, a general feature of UDM models is the possible ap- 
pearance of an effective sound speed, which may become 
significantly different from zero during the evolution of the 
Universe. In general, this corresponds to the appearance of a 
Jeans length (or sound horizon) below which the dark fluid 
does not cluster. Thus, the viability of UDM mode ls strict! 
depends on the value of this effective sound sp eed (|HtjII199; 
IGarriga fc Mukhanovlll999l ; lMukhanovll2005l ) , which has to 
be small enough to allow structure formation (ISandvik et al.l 
|2004 iGiannakis fc Hull2005l ; IBertacca fc Bartoldl2007l ) and 
to reproduce the observed pattern of the CMB temperature 
aniso tropies l|Carturan fc Finellil 120031 ; IBertacca fc Bartold 
120071 ). 

In general, in order for UDM models to have a very 
small speed of sound and a background evolution that fits 
the observations, a severe fine tuning of their parameters 
is necessary. In order to avoid this fine tuning, alternative 
mod els with similar goals have been analysed in the litera- 
ture: iPiattella et al.1 l|2010r i studied in detail the functional 
form of Jeans scale in adiabatic UDM perturbations and in- 
troduced a class of models with a fast transition between an 
early Einstein-de Sitter cold DM-like era and a later ACDM- 
like phase. If the transition is fast enough, these models may 
exhibit satisfactory structure formation and CMB fluctua- 
tions, thus presenting a small Jea ns length even in the case 
of a non-negligible sound speed; iGao et ail l|2009l ) explore 
unification of DM and DE in a theory containing a scalar 
field of non-Lagrangian type, obtained by direct insertion of 
a kinetic term into the energy-momentum tensor. 

Here, we c hoose to invest i gate th e class of UDM mod- 
els studied in Be rtacca et al. (2008), who designed a re- 
construction technique of the Lagrangian, which allows one 
to find models where the effective speed of sound is small 
enough, and the fc-essence scalar fi eld can cluster (see also 
ICamera et al.l [20091 . ICam era 2010). In particular, the au- 
thors require that the Lagrangian of the scalar field is con- 
stant along classical trajectories on cosmological scales, in 
order to obtain a background identical to the background of 
the ACDM model. 

Here, we wish to investigate whether this class of UDM 
models can be scrutinised in realistic scenarios. Specifically, 
we compute the weak lensing signals expected in these mod- 
els as they would be measured by a Euclid-like survey. 

The structure of this paper is as follows. In Section[5]we 
describe the UDM model we use in this work. In Section [3] 
we detail the theory of weak gravitational lensing on the ce- 
lestial sphere, with a particular interest in the cosmic shear 
observable (Section 13. p . In Section [4] we outline the Fisher 
matrix formalism we use to calculate the expected statisti- 
cal errors on cosmological parameters, and with the same 
formalism we compute the expected Bayesian evidence for 
UDM models over the standard ACDM model as a function 
of the sound speed parameter Coo (Section [SJ). In Section [5] 
we present our results, such as the matter power spectrum 
obtained in these UDM models (Section l6.1|) and the corre- 
sponding 3D shear signal (Section 16. 2p ; the parameter esti- 
mations for a Euclid-like survey are presented in Section [6731 
while in Section l"6.4l we use the Bayesian approach to ask the 
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data whether our UDM model is favoured over the ACDM 
model or not. Finally, in Section [Jj conclusions are drawn. 



2 UNIFIED DARK MATTER MODELS 

We consider a UDM model where the Universe is filled with 
a perfect fluid of radiation, baryons and a scalar field <p(t), 
the latter mimicking both DM and DE in form of a cos- 
mological constant. In particular. iBertacca et all (|2008T l. by 
using scalar-field Lagrangians Jzf (X, tp) with a non-canonical 
kinetic term, wherqj 



(i) 



have outlined a technique to reconstruct UDM models such 
that the effective speed of sound is small enough to allow 
the clustering of the scalar field. Specifically, once the ini- 
tial value of the scalar field is fixed, the scalar field La- 
grangian is constant along the classical trajectories, namely 



-A/(8ttG), and the background is identical to the 



background of ACDM. In other words, the energy density of 
the UDM scalar field presents two terms 

/9udm(£) = PDM(i) + PA, (2) 

where pdm behaves like a DM component (pdm oc a -3 ) and 
pA like a cosmological constant component (pa = const.). 
Consequently, Qdm = Pdm (a = l)/Pc and Aa = Pa/ Pc are 
the density parameters of DM and DE today, where p c is the 
present day critical density; hence, the Hubble parameter in 
these UDM models is the same as in ACDM, 



H(z) = H dQ m (l + z 



Qa, 



(3) 



with Ho — 100/ikms -1 Mpc -1 and fi m = Qdm +^6, where 
fit = pt/pc is the baryon density in units of the critical 
density. 

Now we introduce small inhomogeneities of the scalar 
field 5tp{t, x), and in the linear theory of cosmological per- 
turbations and in the Newtonian gauge, the line element is 



ds 2 



-(1 + 2$)dt 2 + a(t)(l + 2*)dx 2 



(4) 



in the case of a spatially flat Universe, as supported by 
CMB measurements (e.g. ISpereel et al.ll2007h . This scalar 
field presents no anisotropic stress, thus = With this 
metric, when the energy density of radiation becomes negli- 
gible, and disregarding also the small baryonic component, 
the evolution of the Fourier modes of the Newtonian po- 
tential <&k(a) are described by ijCarriga fc Mukhanovll 19991 : 
lMukhanovll2005l) 



+ c s k «k 



-«k = o, 



(5) 



where a prime denotes a derivative with respect to the con- 



3 We use units such that c = 1 and signature {—,+,+,+}, where 
Greek indices run over spacetime dimensions, whereas Latin in- 



V PUDM + PUDM ' 
1 



a J 1+ PUDM 
' PUDM 



here, 



2 PUDM,j 

Cs (a) = 



formal time dr = At/a, k — |k| and 

(6) 
(7) 

(8) 

PUDM jX 

is the effective speed of sound, where t x denotes a derivative 

w.r.t. X. 

By following the technique outlined by IBertacca et all 

(2008), it is possible to construct a UDM model in which 
the sound speed is small enough to allow the formation of 
the LSS we see today and is capable of reproducing the ob- 
served pattern of the temperature anisotropies in the CMB 
radiation. We choose a Lagrangian of the form 

&P=pan M (<p,X) = f(<p)g(X)-V{<p), (9) 

with a Born-I nfeld type k i netic term g(X) = 
-VI - 2XM~ 4 l|Born fc Infeldl Il934, l. where M is a 
suitable mass scale. Such a kinetic term can be thought 
as a field theory gen eralisation of the Lagrangia n of a 
relativistic particle tP admanabhan fc Choudhurvl |2002| ; 
lAbramo fc Finellil 120031 ; lAbramo et al.ll2004 ). It was also 
proposed in connection with string theory, since it seems 
to represent a low-energy effective theory of D-branes and 
open strings , and has b een conjectured to play a role in 
cosm ology (|SerJ l2002allbHc3 ; iPadmanabhan fc Choudhurvl 
2002). By using the equation of motion of the scalar field 
ip(t, x) and by imposing that the scalar field Lagrangian is 
constant along the classical trajectories, i.e. pudm = — PA, 
we obtain the following expressions for the potentials 

Acoo cosh(£i/p) 



/(¥>) = 



V(y) = 



sinh(^) [1 + (1 - Coo 2 ) sinh 2 (^)] 



(10) 



l-c 00 2 ) 2 sinh 2 (^) + 2(l 



Coo 2 ) 



1 + (1 



I sinh 2 (£<p) 



(11) 



with £ = V /3A/[4(1 - Coo 2 )!/ 4 ]. Hence, the sound speed 
takes the parametric form 



,(a) = 



SIaco 



Qa + (1 — Coo 2 )S7DMa^ 



(12) 



and it is easy to see that the parameter Coo represents the 
value of the speed of sound when a — > 00. Moreover, when 
a -> 0, c s — > 0. 

In UDM models the fluid which triggers the accelerated 
expansion at late times is also the one which has to clus- 
ter in order to produce the structures we see today. Thus, 
from recombination to the present epoch, the energy den- 
sity of the Universe is dominated by a single dark fluid, and 
therefore the gravitational potential evolution is determined 
by the background and perturbation evolution of this fluid 
alone. As a result, the general trend is that the possible ap- 
pearance of a sound speed significantly different from zero at 
late times corresponds to the appearance of a Jeans length 
ijBertacca fc Bartoldl2007h 



deces label spatial coordinates. 



A j (a) = 



(13) 
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Figure 1. Sound horizon \j(a) for c x = 10" 4 , 10" 3 , 10~ 2 , 1CT 1 
from bottom to top. 



below which the dark fluid does not cluster any more, caus- 
ing a strong evolution in time of the gravitational potential. 
In Fig. [T] we show Aj(a), the sound horizon, for different 
values of Coo. 



3 WEAK LENSING ON THE CELESTIAL 
SPHERE 

In the linear regime, corresponding to the Born approxi- 
mation, where the lensing effects are evaluated o n the null- 
geodesic of the unperturbed ( unlensed) photon (|Hul |200oI ; 
iBartelmann fc Schneider! l200ll ). it is possible to relate the 
weak lensing potential <j> for a given source at a 3D position 
in comoving space x = (%, n) to the Newtonian potential 
$(x) via 



x 



where 



W(x') 



dx 



x-x 



(14) 



(15) 



is the weight function of weak lensing, with n [x( z )] rep- 
resenting the redshift distribution of sources, for which 
Jdxn(x) = 1 holds, and x( z ) being the radial comoving 
distance, such that 

_J_ = dx(£l. (16) 
H{z) dz { ' 

Spin-weighted spherical harmonics and spherical Bessel 
functions are a very natural expansion fo r weak lensing 
observables, such as the potential 4>(x) ( He avens] 120031 ; 
ICastro et al . 2005). Since cosmic shear depends on the New- 
tonian potential, the use of this basis allows one to relate the 
expansion of the shear field to the expansion of the mass den- 
sity field. The properties of the latter depend in a calculable 
way on cosmological parameters, so this opens up the pos- 
sibility of using 3D weak shear to estimate these quantities. 

In the flat-sky approximation, the weak lensing poten- 
tial l[i"4"|l reads 



4>(k,£) 



d x (j)(x)kji (k\) e~ 



(17) 



where I = \£\ is a 2D angular wavenumber, k a radial 
wavenumber and je(kx) a spherical Bessel function of or- 
der I. The covariances of these coefficients define the power 
spectrum of the weak lensing potential via 

(cf>(k, £)(/>* (k',£')} = (2Ti) 2 S D (£-e')C 4 " t '(k,k';e), (18) 

where Sd is the Dirac delta. 

3.1 The 3D shear field 

In this paper we are interested in the information brought 
by the cosmic shear. We now introduce a di stortion tensor 
i|Kaiserlll998l ; IBartelmann fc Schneiderll200lD 



4>,ii{*)= / dx'x'W( X ')<i>^(x,n), 



(19) 



where commas denote derivatives w.r.t. directions perpen- 
dicular to the line of sight. The trace of the distortion tensor 
represents the convergence 



K ( X ) = j (0.1l( x ) +0.22(x)) 

and, defining 71 (x) = i (<^, 



(20) 



L (x) - 0,22(x)) and 72 (x) 



0,i2(x), the linear combination 

7(x) = 71 (x) + i7 2 (x) (21) 

is the differential stretching, or shear. ICastro et all (|2005T ) 
have shown that the complex shear is the second "edth" 
derivative of the weak lensing potential 



7(x) = ^99^ W 



(22) 



where, in Cartesian coordinates {x, y}, 3 = d x + id y . 

We can now express the power spectrum of the 3D cos- 
mic shear as a function of the gravitational potential via 



C 11 {k 1 ,k 2 ;l) = ^ jdkk 2 lf(k 1 ,k)lf(k 2 ,k)P' l '(k,0), (23) 

where P' i (k,z) is the Newtonian potential power spectrum 
and, for a generic field X, we have defined 

I?{ki,k)= [dx^Mw(x)je(k iX ). (24) 



Xk(0) 



4 FISHER MATRIX ANALYSIS 

Cosmological parameters influence the shear in a number of 
ways: the matter power spectrum P s (k,z) is dependent on 
flm, h and the linear amplitude a$. The linear power spec- 
trum is dependent on the growth rate, which also has some 
sensitivity to the parameter of the A-like component equa- 
tion of state wa = Pa/ Pa- It is well know that the speed of 
sound (Eq.[8l is strictly related to wa(z), and it also affects 
the x( z ) relation and hence the angular diameter distance 
sinx [x( 2 )]- These parameters may be estimated from 

the data using likelihood methods. Assuming uniform priors 
for the parameters, the maximum a posteriori probability for 
the parameters is given by the maximum likelihood solution. 
We use a Gaussian likelihood 



21nL 



-TrrinC-C -1 Dl 



(25) 
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where C = ((d — d th )(d — d th ) T ) is the covariance matrix 
and D = (d d th )(d - d th ) T is the data matrix, with d the 
data vector and d th the theoretical mean vector. 

The expected errors on the parameters caiibe esti- 
mated with the Fisher information matrix l|Fisherl [l93i 
Jungman ct al. | Il99rj ; iTegmark et all ll997T l. This has the 
great advantage that different observational strategies can 
be analysed and this can be very valuable for experimental 
design. The Fisher matrix gives the best errors to expect, 
and should be accurate if the likelihood surface near the peak 
is adequately approximated by a multivariate Gaussian. 

The Fisher matrix is the expectation value of the second 
derivative of the InL w.r.t. the parameters {i9 Q }: 



holds; here, c a 



i) 



n + l 



a/3 



d 2 InL 

dd a dd B 



(26) 



and the marginal error on parameter -& a is [(F _1 ) aa ] 5 - If 
the means of the data are fixed, the Fisher matrix can be 
calculated from the covariance matrix and its derivatives 
|Tegmark et al.lll997h by 



"aft 



= — TrfC 1 C t ,C 1 C, 



(27) 



For a square patch of sky, the Fourier transform leads to 
uncorrelated modes, provided the modes are separated by 
27r/0 ra d where ra d is the side of the square in radians, and 
the Fisher matrix is simply the sum of the Fisher matrices 
of each I mode: 



= Q/3 = i£(2£+l)Tr[(c^V Q (dp 



(28) 



where C is the covariance matrix for a given £ mode. 



5 BAYESIAN EVIDENCE 

In this paper we compute parameter forecasts from 3D cos- 
mic shear for UDM models. It is important to notice that 
we are dealing with an alternative model with respect to 
the standard ACDM model; hence, besides determining the 
best-fit value (and the errors) on a set of parameters within 
a model, we can also ask if this particular alternative model 
is preferable to the standard. Model selection is in a sense 
a higher-level question than parameter estimation. While 
in estimating parameters one assumes a theoretical model 
within which one interprets the data, in model selection, 
one wants to know which theoretical framework is preferred 
given the data. Clearly if our alternative model has more 
parameters than the standard one, chi-square analysis will 
not be of any use, because it will always reduce if we add 
more degrees of freedom. From a Bayesian point of view, 
this involves computation of the Bayesian evidence and of 
the Bayes factor B. 

We refer to the two models under examination with 
Mudm and Macdm- We know that, in this context, Macdm 
is simpler than Mudm because it has one fewer parameter, 
i.e. Coo; in the same way, it is also contained in Mudm, be- 
cause, if $ ACDM and i9^/ DM are the parameters of the two 
models (with a = 1, . . . , n and a' = 1, . . . , n + 1), respec- 
tively, then 

K CDM coo} = {^ DM } (29) 



The posterior probability for each model M is given by 
Bayes' theorem 



p(M\d) 



p(d\M)p(M) 



(30) 



The Bayesian evidence is defined as the marginalisation over 
the parameters 



p(d|M) = /d m i9p(d|#,M)p(#|M), 



(31) 



where i9 is the parameter vector, whose length m is n for 
the ACDM model and n + 1 for UDM models. The posterior 
relative probabilities of our two models given the data d and 
with flat priors in their parameters p(M) = const., is then 
(|Heavens et al.ll2007l ; I Heavens! [20091 ') 

p(MACD M |d) _ p(Macdm) 
p(MuDM|d) p(Mudm) 

J d n^ACDM p(d|l? ACDM ;MAcDM)p( ^ACDM |MAcDM) 

X Jd"+ 1 i9 UDM p(dji? UDM ,MuDM)p(i9 UDM |M U DM) ' 

(32) 

If we choose non-committal priors p(Mudm) = 
p(Macdm), the posterior evidence probability reduces to the 
ratio of the evidences, which takes the name of the Bayes 
factor and in the present case reads 



B 



/d"^ ACDM p(d|fl ACDM ,MACDMMfl ACDM |M A cDM) 

J d «+1#UDM p (d|#UD Mi MuDM)p(l? UDM |A/uDM) 



(33) 



Now, let us focus on the priors p(i?|M). If we assume 
flat priors in each parameter, over the range Ai9, then 

p(i9 ACDM |Macdm) = ^ (AtfACDM)-! and 

= fd"/ c ° M P (d|^° M ,M ACDM ) 

J d n+l^UDM p ( d j^UDM )MuDM ) C °°' <> > 

The Bayes factor B still depends on the specific dataset 
d. For future experiments, we do not yet have the data, so we 
compute the expectation value of the Bayes factor, given the 
statistical properties of d. The expectation is computed over 
the distribution of d for the correct model (assumed here to 
be Mudm). To do this, we make two further approximations: 
first we note that B is a ratio, and we approximate (B) by 
the ratio of the expected values, rather than the expectation 
value of the ratio. This should be a good approximation if 
the likelihoods are sharply peaked. 

We also make the Laplace approximation, that the ex- 
pected likelihoods are given by multivariate Gaussians, i.e., 



(35) 



p(d\#,M) = Loe-K*-* )-'*^*-* ^, 

where F a g is the Fisher matrix, given in Eq. 
I Heavens et alj l|2007h have shown that, if we assume that 
the posterior probability densities are small at the bound- 
aries of the prior volume, then we can extend the integration 
to infinity, and the integration over the multivariate Gaus- 
sians can be easily performed. In the present case, this gives 



(B) = 



VdetF UUM L ACDM A Coo 
VdetF AcuM L% DM v^' 



(36) 
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One more subtlety has to be taken into account to 
compute the ratio Lo CDM / Lq DM : if the correct underlying 
model is A/udm, in the incorrect model Macdm the maxi- 
mum of the expected likelihoo d will not, in general, be at the 
correct parameter values (see iHeavens et all |2007| . Fig. 1). 
The n parameters of the ACDM model shift their values 
to compensate the fact that c x is being kept fixed at the 
incorrect fiducial value Coo = 0. With these offsets in the 
maximum likelihood parameters in the ACDM model, the 
Bayes factor takes the form 



= VdetF— ^.. iWjF ^ H 
VdetF ACUM 

where the shifts 5"d a can be computed u nder the assumptio n 
of a multivariate Gaussian distribution l|Tavlor et alj |2007). 
and read 



5&c 



-ACDM \ " 



r UDM r 



(38) 



with G^+j a subset of the UDM Fisher matrix (a vector 
in the present case). 

It is usual to consider the logarithm of the Bayes fac- 
tor, for which the so-called "Jeffreys' scale" gives empirically 
calibrated leve ls of significance for the strength of evidence 
I Jeffreys! 1 196lh . 1 < |ln_Bj < 2.5 is described as "substan- 
tial" evidence in favour of a model, 2.5 < | In B\ < 5 is 
"strong," and |ln_B| > 5 is "decisive." These descriptions 
seem too aggressive: | ln_B| = 1 corresponds to a posterior 
probability for t he less-favoured mode l which is 0.37 of the 
favoured model (|Kass fe Raftery|[l995lV Other authors have 
introduced different terminology fe.g. lTrottall2007l ). 



6 RESULTS AND DISCUSSION 

We use a fiducial cosmology with the following parameters: 
Hubble constant (in units of 100 kms -1 Mpc -1 ) h = 0.71, 
present-day total matter density (in units of critical density) 
Q m = Qum + f^b = 0.3, baryon contribution fif, = 0.045, 
cosmological constant contribution J1a = 0.7, spectral in- 
dex n 3 = 1, linear amplitude (within a sphere of radius 
8 k' 1 Mpc) a 8 = 0.8. 

In Section 16.11 we compute the predicted matter power 
spectrum for UDM models, with a comparison to ACDM. In 
Section [6.21 the 3D shear matrix C 77 (fci, ki\€) is shown. In 
Section f6. 31 we present the parameter forecasts, and in Sec- 
tion [S]3] we show the expected Bayesian evidence for UDM 
models over the ACDM model. 



6.1 The matter power spectrum 

Our class of UDM models allows the value w — — 1 for 
a — > oo. In other words they admit an effective cosmo- 
logical constant energy density at late times. Therefore, in 
order to compare the predictions of our UDM model with 
observational data, w e follow the same prescription used by 
iPiattella etHI (|2010h . where the density contrast of the clus- 
tering fluid is 

8p m PDM<5lJDM + PbSb ,„„, 

d = = , (39) 

Pm Pm 

where St and <5ijdm are the baryon and the scalar field den- 
sity contrasts, respectively, and we emphasise that pdm = 



Pudm — Pa is the only component of the scalar field density 
which clusters. 



6.1.1 Linear regime 

The today matter power spectrum P(k) = P s (k, 2 = 0) 
is the present value of the Fourier transform of the den- 
sity perturbation correlation function. To construct P(k) 
in the ACDM model, we need the growth factor D(z) = 
5(x, z)/5(x, z — 0) on linear scales (i.e. in absence of 
free-streaming) and the transfer function T(k), that de- 
scribes the evolution of perturbations through the epochs 
of horizon crossing and radiaton-matter transition. Here, 
we us e the transfer function suggested by lEisenstein fc Hul 
i|l999l) . which, with an accurate, general fitting formula, 
calculates the power spectrum as a function of the c osmo - 
logical parameters quite efficiently. lEisenstein fc Hul (|l999l ) 
show that baryons are effective at suppressing power on 
small scales compared to DM-only models. Moreover, the 
small-scale limit of this transfer function can be calculated 
analytically as a funct ion of the cosmological parameters 
|Hu fc Eisensteinl Il998l ). Hence, we can write the matter 
power spectrum as 



D(z = 0) 



(40) 



P(k) = 2n 2 S H 2 J T\k) 

here, 5h is a normalisation. 

To obtain P(k) in UDM models, it is useful to remember 
that the class of UDM models we use here is constructed to 
have the same properties of the ACDM model in the early 
Universe; in Eq. which describes the time evolution of 
Fourier modes of the Newtonian potential 4>k(ai), we thus 
set the same initial conditions for both the UDM and the 
ACDM potentials. Gravity is GR, so we can use the Poisson 
equation 



4>k(a) 



2 QmH ° -Ba- 



rn 



which relates $k(a) to the matter power spectrum via 

(5k(a)fc/*(a)) = (2tt) 3 5 d (k - k') P s {k, a). (42) 

Clearly, if we solve Eq. © with c s = 0, we obtain the stan- 
dard ACDM matter power spectrum. 

Fig. [2] shows the matter power spectrum P(k) for 
ACDM and UDM models, for a number of values of Coo. 
By increasing the sound speed, the potential starts to 
decay earlier in tim e, oscillating around zero afterwards 
l|Camera et al.l l2009): at large scales, if Coo is small enough, 
these UDM models reproduce the ACDM model. This fea- 
ture reflects the dependence of the gravitational potential 
on the effective Jeans length A./ (a). It is easy to see that if 
Coo < 10 -3 the perturbations of the UDM reproduce the be- 
haviour of the concordance model within the linear regime 
(the UDM curve for Coo = 10 -3 is virtually on top of the 
ACDM one). Instead, a larger sound speed inhibits structure 
formation earlier in time, thus we observe less power on small 
scales; in this case, the consequence of the oscillatory fea- 
ture of the gravitational potential, due to the non-negligible 
speed of sound, can be clearly seen. 

In principle, the large-scale distribution of galaxies 
could constrain the value of Coo. However, the shape of the 
power spectrum also depends on the normalisation as and 
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Figure 2. Matter power spectra P(k) = P s (k,0) for ACDM 
(solid) and UDM (dot-dashed), with = 10~ 3 , 10 -2 , 10 _1 from 
top to bottom. 



the spectral index n s : therefore, for ei glVGIl Coo as large as 
10~ 2 , an appropriate choice of as and n B can provide a power 
spectrum in agreement with observations, at least on scales 
where non-linear effects are not dominant. In addition, in 
UDM models it is still unclear how the galaxy distribution 
is biased against the gravitational potential of the scalar 
field on small scales. Therefore the large-scale distribution 
of galaxies does not appear to be the best tool to constrain 
this family of UDM models. On the contrary, a weak lensing 
analysis can constrain the matter power spectrum without 
a fine-tuning of either as or the galaxy bias. 



6.1.2 Non-linear regime 

For wavenumbers k > k n \ ~ 0.2 /i Mpc -1 , non-linear con- 
tributions to the evolution of the Newtonian potential (i.e. 
to matter overdensities) become important. In the ACDM 
model, the gravitational potential satisfies Eq. ([5]), but in 
this case c s is the sound speed of the hydrodynamical 
fluid, and therefore can be set equal to zero in the matter- 
domina ted epoch. For c s = , Eq. Q has an analytic so- 
lution jHu fc Eisensteinl Il999l ; liiul 1 20021 ; iMukhanovl 120051 ; 
iBertacca fc Bartold 120071 ) 



*k(a) =A k [l 



H(a) 



da' 



a J H(a') 



(43) 



where the constant of integration is Ak = $k(0)T(fc), with 
T(k) the matter transfer function and <E>k(0) the large-scale 
potential during the radiation-dominated era. 

To perform further calculations on a wider range of 
scales than t hat a llowed by linear theory, we will use the 
ISmith et al.l l|2003h non-linear fitting formulae for P(k) in 
the ACDM model. However, currently there is no linear- 
to-non-linear mapping in UDM models. Nevertheless, as we 
have seen, differences between the ACDM and UDM models 
arise at scales smaller than the sound horizon. With a cross- 
over wavenumber k ~ 1/Aj, if the sound speed is small 
enough to guarantee that A./ is well within the non-linear 
regime we can assume that the non-linear evolution of the 



UDM power spectrum will be similar to the ACDM one. 
A deeper knowledge on this aspect will be the next step of 
the development of UDM models and has to be explored in 
future work. 



6.2 The 3D shear signal 

For a 20,000deg 2 E uclid-like survey (jCimatti et al l 120091 : 
iRefregier et~al]|2010l ). we assume that the source distribu- 
tion over redshifts has the form (jSmail et al.lll994T l 



n(z) oc z e \ z o i 



(44) 



where zq — z m /1.4, and z m — 0.8 is the median redshift 
of the survey. The source number density with photometric 
redshift and shape estimates is 35 per square arcminute. We 
also assume that the photometric redshift errors are Gaus- 
sian, with a dispersion given by a(z) — 0.05(1 + z). 

In order to avoid the high-wav e numb er regime where 
the fitting formulas of ISmith et al.l (|2003l ) may be unreli- 
able, or where baryonic effects might alter the power spec - 
trum (k > 10/tMpc" 1 ; IWhitd [20o3 ; IZhan fc Knox! |2004| ) . 
we do not analyse modes with k > 1.5 Mpc -1 . Note that the 
non-local nature of gravitational lensing does mix modes to 
some degree, but these modes are sufficiently far from the 
uncertain highly no n-linear regime that this is not a concern 
l|Castro et al.ll2005l ). We include angular modes as small as 
each survey will allow, and analyse up to £ max = 5000 (but 
note the wavenumber cut). 

In Fig. [3]we present the 3D shear matrix C 77 (fci, ki\€). 
The first three rows show log 10 C 77 (fci, k 2 ; I) for the ACDM 
model and for a UDM model with = 5 • 10 -4 and = 
5T0 -3 (respectively) in the (fci, fc2)-plane in blue(gray)-scale 
for a number of values of £. In the fourth row we present 
the diagonal elements fc 2 C 77 (fc, k; £) of the 3D shear matrix, 
where the upper (green) curve refers to the smaller speed 
of sound and the lower (green) curve to the greater c^; the 
ACDM (red) curve is virtually on top of the small-Coo UDM 
curve. Finally, in the bottom row we show k 2 times the ratio 
of the diagonal elements C 77 (k,k;£) of UDM models over 
the ACDM model. 

Th e oscillatory feature s of the UDM gravitational po- 
tential (|Camera et alj [20091 ) . whose power spectrum enters 
the shear via Eq. (I23|l . can be clearly seen in the shear sig- 



nal of the UDM model with Coo = 4 ■ 10~ 3 . The bumps in 
the diagonal signal can be easily understood by looking at 
the log 10 C 77 (fci, ki\t) plot, where it is interesting to notice 
how the oscillations take place along any direction, with the 
obvious symmetry along the k\- and fe-axes. Instead, as 
we have noticed in Fig. [2j when the sound speed is small 
enough we do not see any oscillations and the matter power 
spectrum of UDM models is in agreement with ACDM. This 
agreement holds even at non-linear scales k > 0.2 h Mpc -1 . 

Beyond the oscillations, these signals, expected for two 
different values of Coo, show us the effect of the effective 
Jeans length of the gravitational potential. In fact, The 
Newtonian potential in UDM models behaves like ACDM 
at scales much larger than Xj(a) (Eq. I13|) . while at smaller 
scales it starts to decay and oscillate. Hence, at high values 
of £ and k, which correspond to small angular and physical 
scales, respectively, the signal of weak lensing observables, 
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Figure 3. The 3D shear matrix log 10 C 11 (k\, k2',i) for five values of £ in (blue)gray-scale. In the first row we show the ACDM signal, 
while in the second and third rows we present the UDM signal for Coo = 1.0 ■ 10~ 3 and Coo = 5.4 ■ 10 -3 , respectively. The fourth row 
shows the diagonal elements fc 2 C 77 (fc, k; £), and each curve, from top to bottom, refers to the corresponding matrix above. The ACDM 
curve is virtually on top of the small-Coo UDM curve. The fifth row shows the fractional error. 



like cosmic shear, shows the decay of the gravitational po- 
tential. 

Although the UDM signal for = 5 • 10~ 4 appears to 
be in agreement with the ACDM signal (fourth row of Fig. 
[3]), their fractional difference shown in the fifth row is still 
of order unit at k > 1 /iMpc -1 and is not negligible. In fact, 
we will see below in Section 16.41 that this low value of Coo 
still yields a Bayesian evidence which indicates a statistically 
very large difference between this UDM model and ACDM. 

Finally, in Fig. [3] we can also notice that, the higher the 
value of I, the smaller the physical scales are those which 
contribute to the shear signal. This effect is due to the ap- 
proximate Bessel function inequality, I ^ k\, in Eq. (|24[) . As 
the I value increases the diagonal terms of the covariance 
matrix do not become significant until fcxmax ~ l, where 
Xmax = x( 2 max) is the upper limit imposed on the integra- 
tion over the radial comoving distance. 

6.3 Estimation of cosmological parameters 

Once having introduced the method (Section^) and the sur- 
vey design formalism (Section l6.2|) . now we show cosmologi- 
cal parameter forecasts for such a survey and we explore the 



variation in the marginal errors with changes in the sound 
speed parameter Coo. 

By usin g the Fisher matrix analysis outlined in 
I Taylor et al.l ((2007), we calculate predicted Fisher ma- 
trices and parameter constraints for a 20, 000 square- 
degree Euclid-like survey. In all Fisher matrix calcula- 
tions we use a seven-parameter cosmological set {Q m = 
Odm + Qb, £lb, h, Oa, as, n a , Coo} with fiducial values 
{0.3, 0.045, 0.71, 0.7, 0.8, 1.0} for the first six. The Fisher 
matrix is sensitive to Coo, so we compute the evidences at 
twenty fiducial values from 5 • 10 -4 to 5 ■ 10~ 2 . We find 
that the Fisher matrices are unstable for Coo < 10 -3 . This is 
because, when the sound speed is small, the UDM 3D shear 
signal is virtually indistinguishable from that of ACDM, and 
the numerical derivatives w.r.t. Coo thus become unreliable. 

Fig. [4] shows the Fisher matrix elements marginalised 
over all other parameters. In dark blue(gray) we present 
the results for a UDM model with Coo = TO ■ 10~ 3 and 
in light blue(gray) for Coo = 5.4 ■ 10~ 3 . Notice that re- 
sults are shown for universes which are not necessarily 
flat. In non-flat geometries, the spherical Bessel functions 
je(kx) should be replaced by ultraspherical Bessel functions 
(|Heavens et alj|2006h . For the case considered here 
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£ ^> 1 and fc S> (curvatu r e scale) , then (y) — > h (k%) 
|Abbott fc Schaeferlll986l ; IZaldarriaga fc Seliakl |2000| V The 
expansion used is not ideal for curved universes, but it 
should however be an adequate approximation Ri ven cur- 
rent constraints on flatness (e.g. lLarson et alfe lO). 

The Fisher constraints for lensing are large enough 
that for some parameters (as, Qb) the la confidence re- 
gion has an unphysical lower bound. We note that this is 
a symptom of the Fisher matrices Gaussian approximation. 
iTavlor fc Kitchind (|2010l ) address this concern by suggest- 
ing a semi-analytic approach that only assumes Gaussianity 
in particular parameter directions; we leave an implemen- 
tation of this type of parameter error prediction, or a more 
sophisticated likelihood parameter search for future investi- 
gation. 

Before starting the interpretation of such results, it is 
important to underline that what deeply affects the matter 
power spectrum in UDM models, and thus the lensing signal, 
is the presence of an effective Jeans length for the Newto- 
nian potential. Let us focus on Eq. ([5]): we can consider the 
asymptotic solutions, i.e. long wavelength and short wave- 
length perturbations, depending on whether fc <C 1/Aj or 
k 3> respectively. In the former case, the term in 

Eq. ([SJl involving the speed of sound of the scalar field is 
negligible, therefore the solution is formally the same that 
in th e ACDM model (Eq . 1431 , and the Fourier modes $fe (a) 
read (|Bertacca fc Bartoldl2007l) 

H(a) f a da' 



fe<Sl/A 



j (a) oc 



o H{a>) 

instead, in the opposite regime we have 



(fc«l/Aj); (45) 



*fc»i/A J (a) oc 



\/cs(a) 



k I da ■ 



>(«') 



>H(a>) 



(k » 1/Aj). 
(46) 

This means that what enters in the oscillatory dynamics is 
not only Coo , which however plays an important role, but also 
S7dm and Oa, as described in Eq. (j 12p . Therefore, the links 
which connect the expected marginal errors in Fig. U with 
the corresponding fiducial are not quite straightforward. 
Moreover, we find that the Fisher matrix is rather sensitive 
to Coo. The errors we find on the sound speed parameter 
are almost constant, and go from Acoo = 3.0 ■ 10 , for the 
fiducial value = 1.0 • 10 -3 , to Ac^ = 2.6 ■ 10 -5 , when 
Coo = 1.2- 10~ 2 . 

It is already well known that weak lensing can tightly 
constrain the ( Q m , as )-plane, using standard cosmic shea r 
techniques (see iBrown et al.l 120031 ; Semboloni et al.l l2006). 
and 3D weak lensing constrains as in the same way by mea- 
suring the overall normalisation of the matter power spec- 
trum. The expected marginal errors on fi m and as are in 
fact very promising, particularly in the perspective of com- 
bining the cosmic shear dat a with other cosmolog ical observ- 
ables, i.e. CMB or SNela l|Heavens et al.ll2006l ). However, 
the presence of a sound speed can be mimicked in the power 
spectrum, at least in the non-linear regime, by an accurate 
tuning of some para meter values, on top of all as and n s 
(|Camera et al.ll2009l ). This is why the ellipses of those pa- 
rameters get worse for larger values of Coo. 

In UDM models, there is another aspect which is partic- 
ularly interesting to notice: we are able to lift the degener- 
acy between Q m and Qb without using early-Universe data. 




0.002 



0.008 0.01 



Figure 5. Bayes factor — In B for UDM models over the standard 
ACDM model as a function of the sound speed parameter Coo- 



That is because Hdm and Qb enter in the growth of struc- 
tures in two different ways. The expansion history of the 
Universe takes into account only their joint effect, through 
Q m , whereas the speed of sound is determined by Qdm alone. 
In fact we have to keep in mind that in UDM models there 
is a scalar field which mimics both DM and A, but it still 
has proper dynamics different from that of its respective in 
the ACDM model. 



6.4 Model selection 

In Section [5] we showed how the Bayes factor can be used 
to determine which model is favoured by the data. By us- 
ing the Fisher matrix formalism for a Euclid-like survey, we 
compute the Bayes factor B for UDM models over the stan- 
dard ACDM cosmology. We fix flat prior Acoo = 1. 

The large values of — In B derive from the large devia- 
tions 5$ a in Eq. (|38|) which yield an extremely small expo- 
nential. On turn, the deviations <W a are large because, as 
shown in the right-most column of Fig. [4] (i) the ellipsoidal 
confidence regions are narrow, and (ii) they are almost ver- 
tical; in other words, the ACDM parameters that one would 
derive if living in a universe with a non-null Coo would be 
largely biased. 

We conclude that, if UDM is the correct model, there 
would be large evidence for UDM models over ACDM for 
values of Coo > 10 -3 . However, if Coo is so small that the 
UDM peculiar features in the matter power spectrum only 
appear at k 3> 1/tMpc -1 , namely on galactic or smaller 
scales, in principle, we might be unable to distinguish UDM 
from ACDM, unless the non-linear dynamics and/or the ef- 
fects of the baryonic physics on the DM-like dynamics of 
the scalar field are largely different from what we expect in 
ACDM. 



7 CONCLUSIONS 

In this work, we calculate the expected error forecasts for a 
20, 000 square degree survey with median redshift z m — 0.8 
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Figure 4. Expected marginal errors on UDM model cosmological parameters from a 20, 000 deg 2 Euclid-like survey with a median 
redshift z m = 0.8. Ellipses show the la errors for two parameters (68% confidence regions), marginalised over all the other parameters. 
Dark(light) ellipses refer to a UDM model with Coo = 1.0 ■ 10~ 3 (coo = 5.4 • 10~ 3 ). 



such as Euclid <|Cimatti et all 120091 ; iRefregier etTafl l201Ch 
in the framework of unified models of DM and DE (UDM 
models). We focus on those UDM models which are able 
to repr oduce the same Hubble parameter as in the AC DM 
model l|Bertacca fc Bartolol [20071 ; iBertacca et al.ll2008l ). In 
these UDM models, beyond standard matter and radiation, 
there is only one exotic component, a classical scalar field 
with a non-canonical kinetic term in its Lagrangian, that 
during the structure formation behaves like DM, while at 
the present time contributes to the total energy density of 
the Universe like a cosmological constant A. 

In order to avoid the strong integrated Sachs- Wolfe ef- 
fect which typically plagues UDM model s, we follow the 
technique outlined bv lBertacca et all (|2008l ), that allows one 
to construct a UDM model in which the sound speed is small 
enough to let the cosmological structures grow and repro- 
duce the LSS we see today. This can be achieved by pa- 
rameterising the sound speed with its value at late times, 

Coo • 

An effect of the presence of a non-negligible speed of 
sound of the UDM scalar field is the emerging of an effec- 
tive time-dependent Jeans length Xj(a) of the gravitational 
potential. It causes a strong suppression, followed by oscilla- 



tions, of the Fourier modes <l?k(ffl) with fc = |k| > 1/Xj. This 
reflects on the predicted lensing signal, because the latter is 
an integrated effect of the potential wells of the LSS over 
the path that the photons travel from the sources to the 
observer. 

We calculate the 3D shear matrix C 77 (fci, ki\ €) in the 
flat-sky approxi mation for a large nu mber of values of Coo • In 
agreement with lCamera et all (|2009l ). we see that, whilst the 
agreement with the ACDM model is good for small values 
of Coo, when one increases the sound speed parameter, the 
lensing signal appears more suppressed at small scales, and 
moreover the 3D shear matrix does show bumps related to 
the oscillations of the gravitational potential. 

We also compute the Fisher matrix for a Euclid-like 
survey. It has been shown that 3D lensing is a powerful tool 
in co nstrai ning cosmologi c al pa rameters (e.g. ICastro et al.l 
120051 ). and H cavcns et alJ ( 2006) have demonstrated that it 
is particularly useful in unveiling the properties of the dark 
components of the Universe. By using a seven-parameter 
cosmological set {fi m = SIdm + fib, fit, h, Qa, erg, n s , Coo}, 
with one fiducial value for each parameter, except for Coo , for 
which we use twenty values in the range 5T0 -4 . . . 5- 10 -2 , we 
obtain the expected marginal errors. However, the Coo Fisher 
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matrix elements are unstable in the parameter range Coo < 
10 -3 , because the UDM signal is degenerate with respect to 
ACDM. Therefore, we restrict our analysis by considering 
only sound speed fiducial values larger than ~ 10 -3 . We 
get minimal errors that go from Acoo = 3.0 ■ 1CP 5 , for the 
fiducial value Coo = 1.0 ■ 10~ 3 , to Ac^ = 2.6 ■ 10 -5 , when 
Coo = 1.2- 10" 2 . 

In the case of UDM models, 3D lensing is revealed to 
be even more useful for estimating cosmological parameters, 
because since it encodes information from both the geome- 
try and the dynamics of the Universe, it can lift the usual 
degeneracy between the DM and the baryon fractions, Qdm 
and fi;,. This is because in the Hubble parameter, which de- 
termines the background evolution of the geometry of the 
Universe, both Qdm and Qb enter in the usual way, through 
the total matter fraction Q m . On the other side, the speed of 
sound, which affects the structure formation, and thus the 
dynamics of the Universe, is sensitive only on the DM-like 
behaviour of the scalar field, since for baryons c s — holds. 

Finally, we compute the Bayesian expected evidence 
(e.g- lTrottalhoOTT l for UDM models over the ACDM model as 
a function of the sound speed parameter c^. The expected 
evidence clearly shows that the survey data would unques- 
tionably favour UDM models over the standard ACDM 
model, if its sound speed parameter exceed ~ 10 -3 . 
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